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ABSTRACT 

We examine the performance of four different methods which are used to measure 
mass segregation in star-forming regions: the radial variation of the mass function 
A^mf; the minimum spanning tree-based Amsr method; the local surface density 
Sldr method; and the Ogsr technique, which isolates groups of stars and determines 
whether the most massive star in each group is more centrally concentrated than the 
average star. All four methods have been proposed in the literature as techniques for 
quantifying mass segregation, yet they routinely produce contradictory results as they 
do not all measure the same thing. We apply each method to synthetic star-forming 
regions to determine when and why they have shortcomings. When a star-forming 
region is smooth and centrally concentrated, all four methods correctly identify mass 
segregation when it is present. However, if the region is spatially substructured, the 
Hgsr method fails because it arbitrarily defines groups in the hierarchical distribu¬ 
tion, and usually discards positional information for many of the most massive stars 
in the region. We also show that the Amsr and Sldr methods can sometimes pro¬ 
duce apparently contradictory results, because they use different definitions of mass 
segregation. We conclude that only Amsr measures mass segregation in the classical 
sense (without the need for defining the centre of the region), although Sldr does 
place limits on the amount of previous dynamical evolution in a star-forming region. 
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1 INTRODUCTION 

Most young (<10Myr) stars are observed in the company 
of others, either in clusters dbada fc Lada 2 00^ . groups 
dPorras et al.l l2003l l , associations ( Blaauwl Il964l l or in re- 

? ions which are over- dense with respect to the Galactic field 
Bressert et ahlkoiol l. Attempts have been made at distin¬ 
guishing isolated versus cl ustered star-formation, or clus¬ 
ters versus associations (e .g. IC artwright Sz Whitworth|[2004l : 
iGieles fc Porteeies Zwara I 2 OIII I. However, it appears that 
(at least on local scales ) there is no fundam ental scale unit 
for star formation (e.g. iBressert et al. l l201Cll ). and that clus¬ 
ters may just b e the dense tail o f the distribution of Galactic 
star formation (lKruiisserJl2012l ). 

Given that star formation on pc-scales often seems to 
produce complex, hierarchical distributions, it is crucial to 
be able to examine the spatial distributions of stars in a 
quantitative, statistical way. This allows us to extract infor¬ 
mation on star formation and compare and contrast different 
regions (and simulations). One such indicator is the relative 
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positions of the most massive stars with respect to the aver¬ 
age stars in a region - often referred to as ‘mass segregation’ 
when the most massive stars are more centrally concentrated 
than average. 

In recent years several alternative methods for measur¬ 
ing mass segregation have been proposed in the literature. 
In this paper we critically ass ess several of the ne w methods, 
including the Am sr method llA llison et al.||2009l). the E — m 
or Eldr method i Maschberger fc Glarkell201ll ) and the tech¬ 
nique of dividing a star-forming region into groups to deter¬ 
mine the relative dis tance of the most massive star to th e 
centre of each group llKirk &: Mverj[201ll : lKirk et al.ll2014l ). 
which we call Hgsr- A similar study to compa re mass segre- 
jation methods was previously carried out bv lOlczak et al.l 
2 OIII) : however, the Eldr and Hgsr methods were not yet 
prominent in the literature and therefore not included in 
that study. A full appraisal of these methods - and Amsr “ 
using synthetic data is therefore required in order to make 
sense of the growing literature in this area. 

The paper is organised as follows. In Section [5] we pro¬ 
vide a brief summary of the different methods used to define 
mass segregation, and discuss the methods we will use in this 
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paper in detail. In Section Owe test these methods on syn¬ 
thetic data, before discussing our results in Section [4] and 
concluding in Section O 


2 WHAT IS MASS SEGREGATION? 


‘Mass segregation’ is a phrase often used in relation to star 
clnsters and regions, but often it is not defined. 

The classical definition of mass segregation is based on 
the behaviour of dynamically old bound systems (i.e. relaxed 
and virialised star clusters). The process of two-body relax¬ 
ation redistributes energy between stars and they approach 
energy equipartition whereby all stars have the same mean 
kinetic energy - this means that more massive stars will have 
a velocity dispersion that is lower. Because the velocity dis¬ 
persion of the more massive stars is lower, they will tend to 
be co ncentrated towards the centre of the clnster dSpitzeil 

The timescale, tMs{M), on which stars of mass M will 
dynamically mass segregate depends on the mean mass of 
stars (m) and the two-body relaxation time, treiax llSpitzeil 

Il969lf 


tmsiM) 


M 

(m> 


^relax • 


( 1 ) 


Therefore in old star clusters (especially globular clusters) 
we expect to see mass segregation down to low masses and 
that mass segregation is explicable entirely by dynamics. 

In these old clusters, the two-body dynamics that 
may have caused any mass segregation will also remove 
any primordial substructure in the spatial distribution and 
we would expect (and observe) the cluster to have a 
s mooth, centrally co ncentrated sp atial distribution, such as 
a lPlummerl dlQllIl or iKind (Il966l l profile. 

In this case, the clusters have a well-defined radial pro¬ 
file where one can quantify mass segregation by taking dif- 
ferent mass bins and comparing the density profiles (e.g. 
iHillenbrandl 1 19971 : IPinfield et al.l1l998l i , or examining varia¬ 
tions in the slope of the mass function (or luminosity func- 


tion) with distance from the cluster centre ( 

CarDenter et alJ 

1 19971: Ide Griis et al.ll2002l: IGouliermis et all 

2004|). A related 


method is to quantify the variation of the ‘Spitzer radius’ - 
the rms distance of star s in a cluster around t he centre of 
mass - with luminosity (iGouliermis et ai1l2009li . 

The motivation in attempting to observe mass segrega¬ 
tion in young clusters or regions is that it might not have a 
dynamical origin. If we observe mass segregation in a region 
that is so young that two-body encounters cannot have mass 
segregated the starfl, then the mass segregation must be set 
by some aspect of the star formation process, and is often 
labelled ‘primordial mass segregation’. 

Primordial mass segregation h as been found in some 
simu l ations of star form ation (e.g. iMaschberger fc Clarke 


20121: Parker et al. 


20 111 : lM^ers_et_^ 201^ but not in others ( Girichidis et al 


2015l b and we also note that any ob¬ 


served signature may be a combination of primordia l and 
dynamical mass segregation (iMoeckel fc Bonnellll^03 1. and 


^ This is complicated by the fact that it is not the current re¬ 
laxation time of the region that is important, rath er it is how 
dynamically old the region is Isee lAllison et aLll201Clh . 


to complicate matters further mass segregation can be intro¬ 
duced very rapidly dynamically thro ugh violent relaxat ion 
rather than two-body relaxation 1see I Allison et al]l2010ll . 


Observationally, mass segregation has been searched 
for in clu sters and star f orming regions for many 


years (e.g. ISagar et all Il988l: Ide March! & Parescd 

1996: 

Hillenbrand & HartmanrJI 19981: Raboud & Mermilliod 

1998: 

de Griis et al.l 20021: Gouliermis et al. 2004: Sabbi et al. 

20081: IGouliermis et al.||2009|. and manv morel, but in the 


past few years a number of new statistical methods have 
been developed for finding mass segregation, and the pur¬ 
pose of this paper is to examine what exactly they are search¬ 
ing for, and what problems they might have. 

Here we briefly describe the most commonly used meth¬ 
ods and their assumptions. In the next section we describe 
exactly how we implement each method in detail. 

It should be noted that all methods suffer from a poten¬ 
tially very serious problem. All methods examine the relative 
distributions of high- and low-mass stars. Therefore the dis¬ 
tribution of low-mass stars must be known, however these 
are very faint compared to the high-mass stars in young re¬ 
gions/clusters. And so the location in which the observer 
is biased against observing low-mass stars is near luminous 
high-mass stars which is exactly where one needs to know if 
low-mass stars are present or not. In this paper we use fake 
data in which we have the advantage of knowing exactly 
where every star is and what its mass is, but real obser¬ 
vations do not hav e this advantage (see IStolte et al.l l2005l : 
lAscenso et al ]|2009l. for a detailed discussion of observational 
selection effects and biases). 


2.1 Radial mass functions, Mmf 


Until recently, the most commonly used way of determining 
mass segregation was to compare the radial distribution of 
the mass function, or the radial distribution of stars of dif- 


j.lSagar et al.l 

1984 

1: ISabbi et al.l 

2008|: 


Ghavarn'a et al. 20101 1. 


In this definition, if a cluster is mass segregated the most 
massive stars are preferentially towards the centre, and low- 
mass stars preferentially in the outskirts. Therefore the slope 
of the mass function should be flatter in the centre than the 
outskirts, and/or the low-mass stars should have a much 
broader radial distribution. 

This method suffers three significant drawbacks. One is 
that it requires a ‘centre’ to be defined. In relaxed, virialised 
star clusters there is a centre from which this can be mea¬ 
sured, but in substructured and messy young star forming 
regions it is unclear if defining a ‘centre’ means anything at 
all, and even if it did, it is unclear how one would practically 
do this. One can define a geometric centre based on the po¬ 
sitions of stars, but in some morphologies (as we shall see) 
this ‘centre’ is not near to, nor is the central location, of the 
stars. 

Second, the definition of radial bins is non-trivial and 
often arbitrary, which adds a further level of complexity to 
the interpretation of any signal. 

It also suffers from poor statistics, in that it is unclear 
how to compare different clusters or provide any quantitative 
information on the mass segregation. 

In practice, we compare the cumulative distribution of 
the radii of the ten most massive stars with the cumulative 
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Comparative mass segregation 3 


distribution of the radii of all the stars in the distribution. 
We set the ‘centre’ of the region to be the (known to us) 
centre of mass of the region. We gauge the significance of 
any difference using a Kolmogorov-Smirnov (KS) test, and 
rather generously use a p-value of <0.1 as our significance 
threshold. 


2.2 The AMSR-parameter 


In order to try and avoid the problems of defining a centre 
and to produce a quantitative statistic to aid comparisons 
Allison et al. (2009) introduced the AMSR-paranieter. 

The AMSR-parameter examines the relative (2D) spatial 
distributions of the N most massive stars relative to each- 
other with the spatial distributions of N random stars. This 
is done by finding the edge length of the minimum spanning 
tree (MST) that connects the N most massive stars with 
many MSTs of N random stars. 

This has the great advantage that it does not require a 
centre or any special position to be dehned. It also produces 
a number with associated error that states how much longer 
or shorter the massive star MST is compared to random 
MSTs, and how likely it is that this length could be drawn 
at random from the random MSTs (i.e. how likely is this to 
occur by random chance). 

The AMSR-parameter measures a very similar ‘mass seg¬ 
regation’ to the classical definition: it examines if the mas¬ 
sive stars are closer to each-other than one would expect 
by random chance. If the massive stars are closer to one 
another, then the star-forming region is said to be mass seg¬ 
regated. 

It practice, we take the ratio of the average (mean) ran¬ 
dom MST length to the subset MST length, a quantitative 
measure of the degree of mass segregation (normal or in¬ 
verse) can be obtained. We hrst determine the subset MST 
length, /subset. We then determine the average length of sets 
of A’mst random stars each time, (/average). There is a dis¬ 
persion associated with the average length of random MSTs, 
which is roughly Gaussian and can be quantified as the stan¬ 
dard deviation of the lengths (/average) ± craverage- However, 
we conservatively estimate the lower (upper) uncertainty as 
the MST length which lies 1/6 (5/6) of the way through an 
ordered list of all the random lengths (corresponding to a 66 
per cent deviation from the median value, (/average))- This 
determination prevents a single outlying object from heav¬ 
ily influencing the uncertainty. We can now dehne the ‘mass 
segregation ratio’ (Amsr) as the ratio between the average 
random MST pathlength and that of a chosen subset, or 
mass range of objects: 


Amsr = 


// \ "tUg 

Xt-average / 

/subset 


( 2 ) 


A Amsr of 1 shows that the stars in the chosen subset are 
distributed in the same way as all the other stars, whereas 
Amsr > 1 indicates mass segregation and Amsr < 1 indi¬ 
cates inverse mass segregation, i.e. the chosen subset is more 
widely distributed than the other stars. 

Note that a slightly different formulation of Amsr, 
which uses the geometric mean to define the uncertainties, 
is available (lOlczak et al.ll2011^. but for the purpos es of this 
paper the original method by lAllison et al.l (l2009ll is suffi¬ 
cient. 


2.3 The local density ratio, Eldr 

iMaschberger fc Clarkd l^2011^ introduced another measure 
of mass segregation - the local density ratio. For this the 
local surface density of every star E is found, and the average 
local surface density of the N most massive stars is compared 
to the average local surface density of all stars to obtain 
the ‘local surface d ensity ratio’, Eldr (iKupoer et al.|[201ll : 
IParker et al.l[20l3l . This is able to determine if the most 
massive stars are in regions of significantly higher (or lower) 
local surface density than would be expected by random 
chance. If the most massive stars are in areas of higher local 
density than the region average, the region is said to be mass 
segregated. 

As with Amsr, Eldr makes no assumptions about there 
being a centre. However it is possible that the same numeri¬ 
cal value of Eldr can be statistically signihcant sometimes, 
but not other times, which we quantify by means of a KS 
test. 

It is very important to note that Amsr and Eldr mea¬ 
sure different versions of ‘mass segregation’. Amsr deter¬ 
mines if the massive stars are closer to each other than one 
would expect by random chance, Eldr determines if the 
massive stars are in regions of higher surface density than 
one would expect by random chance. 

It would be quite possible for Eldr to hud ‘mass segre¬ 
gation’, but for Amsr not to. This would occur for example 
if the massive stars were widely distribu ted, but each had a 
local overdensity of low-mass stars (see IParker et al.l l2014l . 
for examples). 

In practice, we calculate the local stellar surface den - 
sity following the prescription of ICasertano fc Hud lll985h . 
modihed to account for the analysis in projection. For an 
individual star the local stellar surface density is given by 


E = 



(3) 


where is the distance to the n*** nearest neighbouring 
stai0. 

iKiipper et al.l (l201ll l and IParker et all (|2014h took the 
ratio of the median surface density of a chosen subset (in 
this paper we will use the 10 most massive stars) to the 
median for all stars in a region to dehne the local surface 
density ratio, Eldr: 


Eldr 


Esubset 


(4) 


The Eldr ratio is then quoted with the p-value from the 
KS test to gauge the signihcance of any deviation from the 
median for all stars, again with a p-value < 0.1 used as the 
boundary between the difference being signihcant or not. 


2.4 Group segregation ratio, Hgsr 

A further, alternative method fo r quantifying mass s egre¬ 
gation was recent ly suggested bv iKirk fc Mver3 ll201lll and 
iKirk et ^ (l2014h . In this method stellar groups are identi- 
hed and ‘isolated’ from the total distribution. Each of these 

^ Note that this n does not need to have the same value as the 
N used to define the subset of choice, although we adopt n = 
10, A = 10 throughout this work. 
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groups is then examined to see if the most massive star it 
contains is closer to the centre of the group than the median 
distance of all stars in the group. If the most massive star 
is closer to the centre than the average star in the major¬ 
ity of the groups, the star-forming region is said to be mass 
segregated. 

The ‘mass segregation’ searched for in this method is 
different again from both Amsr and Eldr. flcsR divides 
the region into groups and then examines each group for 
evidence of internal mass segregation. In this process many 
stars in the region (possibly including some of the most mas¬ 
sive) can be excluded if they do not belong in a group. 

First, an MST is constructed for the entire region and 
a cumulative distribution of all MST branch lengths is then 
made. Two power-law slopes are then fitted to the short¬ 
est lengths, and the longest lengths, and the intersection of 
these slopes defines the boundary of subclustering, dbreak 
llGutermuth et alll2009ll . 

The links in the full MST which exceed dbreak are re¬ 
moved, resulting in the star-forming region being divided 
into groups. If the position of the most massive star in the 
group rmm is closer to the central position than the median 
value for all stars, Tmed, the group, or subcluster has an off¬ 
set ratio (rmm/rined) less than unity and is said to be mass 
segregated. 

In this paper, we define a ‘group segregation ratio’, 
f^GSR, as 

f^GSR = (5) 

iVgrp 

where A^grp are the number of groups, and A^aeg is the number 
of these groups that have an offset ratio less than unity. If a 
star-forming region has Hgsr = 1, then all individual groups 
defined by dbreak are mass segregated. 

It should be noted that, just by random chance we 
would expect Hgsr ~ 0.5, as half of the time the most mas¬ 
sive star would be in the inner 50 per cent of stars. The 
significance of Hgsr ~ 0.5 is affected by Poisson noise; for 
example an Hgsr = 0.8 is not significant if it is 8 out of 10 
subgroups. 


3 FINDING MASS SEGREGATION IN 

SIMULATED REGIONS 

All of these four methods for finding ‘mass segregation’ will 
find classical mass segregation in a relaxed, spatially smooth, 
spherical, bound cluster. If the most massive stars are close 
together in the centre of a spherical cluster then: (A) Mmf 
will show a different mass function in the inner regions. (B) 
Amsr will show that the massive stars are concentrated to¬ 
gether. (C) Eldr will find that the most massive stars are 
in the regions of highest surface density. (D) Hgsr will find 
that the most massive star is towards the centre of a single 
group (in this situation, the cluster itself is the group). 

In such a situation we would advise the reader to use 
Amsr to look for mass segregation as it gives a single quanti¬ 
tative value for the degree of mass segregation with an asso¬ 
ciated error, and can easily determine the stellar mass down 


to which mass mass seg r egation is present (e .g. Allison_et_^ 


I2OO9I : ISana et ^ I2OIC ; iBeccari et al.l I2OI2I : iDelgado et al 


201^ E£_e£^ 2OI3I : IPang et al.l I2OI3I : iRivilla et al.l l2014l : 
Wright et al. boi4i. 


What we will examine in this paper is the analysis of 
complex, substructured regions and the search for ‘mass seg¬ 
regation’ within them, according to the definition presented 
in each method. 

In this Section we run our four methods for quanti¬ 
fying mass segregation on a synthetic dataset containing 
N = 300 stars to match the small-A statistics of many 
young regions. We distribute the st ars in a fractal distribu¬ 
tion a cc ording to the p rescri ption in Goodwin fc WhitworthI 
||2004) , lAllison et al.l (l2010l ~) and I Parker et al.l 1 20l4) . "VW 
refer the reader to those papers for a detailed description 
of how the fractal is set up, but we briefly summarise the 
method here. The fractal is built by creating a cube con¬ 
taining ‘parents’, which spawn a number of ‘children’ de¬ 
pending on the desired fractal dimension, D. The amount 
of substructure is then set by the number of children that 
are allowed to mature (the lower the fractal dimension, the 
fewer children mature and the cube has more substructure). 

We choose a fractal distribution because the Amsr, 
Eldr and Hgsr methods were developed specifically to 
be used on substructured, or hierarchical spatial distribu¬ 
tions of stars in star-forming regions and clusters. All three 
negate the requirement of designating a ‘centre’ (although 
the flGSR method requires the definition of group centres 
- as detailed in Section 12.411. Star-formation may result in 
a truly fractal distribution l|Elmegreen fc ElmegreenI [ 200 II : 
iGartwright fc Whitworthl2004^ . but star forming regions are 
unlikely to retain their primordial spatial d istribution due 
to dynamical evolution llParker et al.l 120141 . As a default, 
we choose D = 2.0 and assign the fractal a radius of 5 pc. 

We note that observed sta r-forming regions display 
a ra n ge of fractal dimen s ions (lCMliWright_£^(Wlhityrorth| 
2004 ^hmeia e t al.l l2008l . I2009I : ISanchez fc Alfarol 20091 : 


Gouliermis et al.ll2014^ . Tt is possible that all regions may 


form with the same (low) fractal dimension, and ob¬ 
served diff erences are due to differing amounts of dynamica l 
evoluti on llGoodwin fc WhitworthI l2004l : IParker et al.ll20i3 : 
[Parked l2014h . although simulations of star-formation also 
produce a range of values before significant dynamical evo- 
lution takes place (iGirichidis et al.ll2012l : [Dale et ^ l2012l . 
l2013li . For the purposes of the numerical tests presented 
here, regions with D = 2.0 adequately highlight the differ¬ 
ences between the various defi nitions of ma s s segr egation. 

We draw masses from the iMaschberged ||2013|) formula¬ 
tion of the initial mass function (IMF): 


p(m) oc ( — 

n 


1-f ( — 

n 


( 6 ) 


Eq. [6] essentially combin es the log - norm al app roximation 
for the I MF d erived by IChabried (l2003l . l2005h with the 
ISaloeted lll955ll power-law slope for stars with mass >1 Mq. 
Here, /r = 0.2 M© is the average stellar mass, a = 2.3 
is the Salpeter power-law exponent for higher mass stars, 
and P = 1.4 is the power-law exponent to describe the 
slope of the IMF for low-m ass objects (which also deviate s 
from the log-normal form; iBastian. Govev fc Meveill2010l ~l. 
Finally, we sample from this IMF within the mass range 
miow = 0.01 M 0 to ruup = 50 Mq. 

We note that in this paper, the choice of the mass distri- 
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bution is unimportant, because we are comparing the posi¬ 
tions of the 10 most massive stars to the positions of all stars 
in a region. In this case, we only require the massive stars 
to have masses above those of the remaining 290 objects. 

However, stochastic sampling of the IMF for low-A'’ re¬ 
gions can res ult in very different mass distributions between 
realisations (IParker fc GoodwinI 1200^ 1. and if we were to 
study the subsequent dynamical evolution, differences in the 
relative masses of individual stars could affect the measured 
amount of mass segregation. 

We created 20 realisations of the fractal distribution 
(and mass distribution). However, in the following we con¬ 
centrate on just one typical realisation of the spatial dis¬ 
tribution, and mass distribution of stars. The realisation is 
‘typical’ in the sense that it nicely highlights the advantages, 
and disadvantages of each of the methods. 

For each realisation we assign the stellar masses in one 
of three ways. First the stellar masses are randomly assigned, 
so no mass segregation - whatever the definition - should be 
detected ('Section l3.1ll . We then change the positions of the 
ten most massive stars so that they are either more centrally 
concentrated fSection 13.21) or in the locations of the highest 
stellar surface density fSection 13.31) . 

We perform all the subsequent analysis in 2 dimensions 
in order to mimic the information available to observers. 

It should be noted that all methods make a hidden as¬ 
sumption that the 2 dimensional distributions are a good 
representation of the true 3 dimensional structure (in the 
sense that they retain the same information on spatial dis¬ 
tributions). It is unclear if this is really the case, and we will 
return to this difficult question in a later paper. For now we 
will proceed under the assumption that the 2 dimensional 
distributions do retain the important information present in 
the true 3 dimensional structure. 


3.1 Random distributions of stellar masses 

In this section we use our four techniques for measuring mass 
segregation on a random distribution of stars, as shown in 
Fig.Hl The locations of the ten most massive stars are shown 
by the large red dots. 


3.1.1 Radial mass functions, AImf 


In Fig. 2(a) we show the cumulative distribution of the ra¬ 


dial distances from the centre of the fractal for all stars (the 
solid black line) and for the ten most massive stars (the red 
dashed line). Within 3pc of the centre, the two distribu¬ 
tions are overlaid. However, there are no massive stars at 
radii greater than 3.2 pc and the cumulative distributions 
differ beyond this radius. However, a KS test on the two 
distributions returns a p-value of 0.51 that the two popula¬ 
tions share the same parent distribution, i.e. the difference 
is not significant in that it is higher than our threshold of 
p = 0.1. 

It is important to note that the centre of the distribution 
is known to us to be at {0, 0} pc. When confronted with a 
distribution similar to that shown in Fig. [T] an observer 
would use the average position of all the stars to define a 
centre. In the example shown here, this average position is 
almost identical to the centre of mass. 


4 
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(J 

Q_ 
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Figure 1. A fractal distribution (D = 2.0) with stars randomly 
drawn from an initial mass function and placed randomly in the 
spatial distribution. The ten most massive stars are shown by the 
larger (red) points. 
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3.1.2 Mass segregation ratio, Amsr 


In Fig. |2(b)] we show Amsr as a function of the A^mst stars 
for the fractal region in Fig [T] Amsr = 1 (consistent with 
there being no mass segregation) is shown by the horizon¬ 
tal red dashed line. The Amsr technique shows that the 10 
most massive stars are slightly more centrally concentrated 
than the average stars, with Amsr = ITIq 4 stars with 
m > 2.05. The 20 most massive stars are also slightly more 
centrally concentrated than the average stars. 

This positive signal of mass segregation is likely due to 
the same spatial feature in Fig. [T] that shows an apparent 
difference in the radial mass functions, namely that none of 
the most massive stars are more than 3 pc from the centre 


(Fig. 2(a)). This is a 2-a difference from unity, and so would 
be expected roughly l-in-20 times. By ‘fluke’ this is the only 
random realisation that shows a 2-a signature of mass seg¬ 
regation and emphasises the need to avoid over-interpreting 
a single 2-a result. 


3.1.3 Local density ratio, Eldr 


In Fig. 2(c) we plot the local stellar surface density, E, 


against individual stellar mass m in Fig. 2(c) The me¬ 


dian stellar surface density for the entire distribution is 
13.1 stars pc“^ (the blue dashed line) whereas the median 
stellar surface density for the ten most massive stars is 
13.2starspc“^ (the solid red line). These values are obvi¬ 
ously very similar, and a KS test returns a p-value of 0.3 
that they share the same parent distribution - i.e. this is 
not a significant difference compared to our threshold value 
of p = 0.1. Eldr = 1.0 and we therefore conclude that the 
most massive stars are not mass segregated according to this 
definition. 
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Figure 2. Three separate measures of mass segregation for the stellar distribution shown in Fig.[T] In panel (a) we show the cumulative 
distribution of the distance from the centre for the ten most massive stars (the red dashed line) and the cumulative distribution for all 
stars (the solid black line) - the mass function comparison, A4mf- In panel (b) we show the Amsr mass segregation ratio as a function of 
the A^mst stars used in the subset (the lowest mass star, m^, for various A^mst values is shown along the top horizontal axis). Amsr = 1 
(i.e. no preferred spatial distribution) is shown by the solid horizontal red dashed line. In panel (c) we show local stellar surface density 
versus stellar mass (the S — m plot). The median stellar surface density for the ten most massive stars is shown by the righthand solid 
(red) horizontal line and the median surface density for all of the stars is shown by the blue horizontal dashed line. 


3 . 1.4 Group segregation ratio, flcsR 

In Fig. |3] we show the determination of the Hgsr group 
segregation ratio. We start by constructing an MST of the 
entire region, as shown in Fig. 3(a) The cumulative dis¬ 


tribution of the branches in the entire MST is shown in 


Fig. 3(b) and the two power law fits to the close branches 


and th e long branches ar e show n by the solid re d lines . Fol- 
lowing iG uterm uth et al] (l2009l ): iKirk &: MversI ll201ll) and 
iKirk et al.l 1 ~14h we take the intersection of these lines as 
the MST dbreak = 0.28 pc, the boundary between ‘clus¬ 
tered’ and ‘diffusely’ distributed stars. All MST branches 


with length < dbreak are retained (Fig. 3(c) I which defines 
groups within the distributior0. 


In Fig. 3(c)| we show the groups (with N > 2) selected 
by this method as the stars still connected by the red MST 
links. The most massive star within each of the groups is 
marked by a red triangle. The ten most massive stars in the 
entire region are shown by the large open blue circles. 

It might be considered that dbreak = 0.28 pc has a physi¬ 
cal importance - it is the apparent break between structures. 
However, in this simulation this distance has no physical 
significance, it is just a projected 5 pc radius D = 2.0 box 
fractal distribution, which by design is hierarchical and self¬ 
similar, with no special spatial scale (apart from the radius 
itself). Examination of Fig. |3] shows that there is nothing 
‘special’ about this distance. Furthermore, small changes to 
dbreak Can drastically affect the number of groups which are 
identified. For example, if we choose dbreak = 0.25 pc (the 


point at which the cumulative distribution in Fig. 3(b) de¬ 


viates from the steep power-law), we identify 18 groups (in¬ 
stead of the 16 identified using dbreak = 0.28 pc). If we 
choose dbreak = 0.5 pc (where the shallow power-law slope 
deviates from the tail of the distribution), we identify only 
6 groups. 


® Other methods to define groups/clusters in cr owded fiel d s may 
have advantages over the MST technique - see ISchmei3 i 
for a review. 


Examination of Fig. [T] shows to the eye perhaps five 
groups, the most significant being to the bottom right. 


Fig. 3(c) shows that Hgsr has identified many more groups 


than this. In particular, the stars to the bottom left (around 
{—3, —1} pc) have been split into three groups. And the sig¬ 
nificant distribution of stars at the bottom centre/right have 
been split into several groups, but some stars (including one 
of the most massive in the region) have not been included 
in any group. 

The identification of groups is crucial to the Hgsr 
method, but it is unclear from Fig. |3(c)| that the selected 
groups are ‘real’ in any sense. 

There are three other significant issues with the Hgsr 


method that are immediately apparent from Fig. 3(c) 


Firstly, the most massive star in a group may not be one 
of the most massive stars in the region. All of the groups to 
the upper left have a locally most massive star that is not 
one of the most massive stars in the region as a whole. 

Secondly, a group may contain more than one of the 
truly most massive stars in a region (e.g. three of the larger 
groups to the bottom right) in which case only the most 
massive of these is considered and the other (truly massive 
for the region) stars are discarded. 

Thirdly, if a truly massive star is not part of a group 
(surely an interesting phenonena) then it is discarded from 
the analysis entirely (e.g. the large open blue circle at the 
bottom centre). 

But taking the Hgsr method to its conclusion we show 
the mass ratio (most massive star in the group to median 
star in the group) versus offset ratio (position of most mas¬ 
sive star to median group position) in for all groups with 
N > 2 stars in Fig. 3(d) The five groups with ^ 10 stars 
are shown by the red points, and three of them have an offset 
ratio less than unity, i .e. they are mass seg rega ted accord¬ 
ing to the definition in iKirk fc Mverj (1201 ll ) and iKirk et al.l 
(|2014) . We define the group segregation ratio as the num¬ 
ber of groups with an offset ratio less than or equal to 
unity divided by the total number of groups. For all groups 
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(d) Mass ratio versus position offset ratio for groups 


Figure 3. Mass segregation as defined by the Qgsr method. In panel (a) an MST of the full spatial distribution is shown, and the 
cumulative distribution of all of the branch lengths is shown in panel (b). The two power law slopes used to fit the data are shown by 
the red lines. The intersection of these slopes gives the critical MST length, dbreak^ m panel (c) we show the groups identified using 
this length. In the groups in which there are 3 or more stars the most massive star in the group is shown by the solid red triangle. The 
positions of the ten most massive stars in the full distribution are shown by the large open blue circles (these correspond to the filled 
red circles in Fig. 0. In panel (d) we show the mass ratio of the most massive star in each group to the group median mass versus the 
ratio of the position of the most massive star to the median position of the group. Groups with ten or more stars are shown by the red 
asterisks. 


with N > ^r^GSR = 0.59 and for groups with N ^ 10, 
f^GSR = 0.6Qj. In a truly random distribution, f^GSR = 0.5. 

It is worth noting another problem here, in that the 
f^GSR method needs to define a ‘centre’ of each group from 
which to measure distances. Therefore there is an implicit 

Note that iKirk &; Mverj ll201lh and iKirk et al.l |2o3) generally 
only present statistics for groups containing 10 or more stars, and 
we will also draw conclusions based only on these ‘large’ groups. 


assumption of spherical symmetry which examination of 
Fig. |3(c)] shows not to be the case in most groups. 

There are two groups with ^ 40 stars, which is 
enough points to run the Amsr and Eldr methods on these 
groups in isolation. One of these groups has an offset ratio 
of less than unity (i.e. it is mass segregated according to this 
method) - however, neither Amsr nor Eldr find mass seg¬ 
regation in this group. The second group has an offset ratio 
of greater than unity (i.e. it is inversely mass segregated ac- 
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Figure 4. As Fig. [T] a fractal distribution of stellar masses ran¬ 
domly drawn from an initial mass function. However, in this case 
we have swapped the locations of the ten most massive stars 
(shown by the larger red points) with the ten most central stars. 



cording to this method). Again, both Amsr and Eldr are 
consistent with no mass segregation (normal or inverse). 


3.2 Massive stars centrally concentrated 

We now swap the positions of the 10 most massive stars 
with the positions of the 10 stars closest to the centre of the 
fractal distribution, as shown by the red points in Fig. U) 
This is clearly a rather artificial distribution of the most 
massive stars, but it is one that most closely matches the 
‘classical’ definition of mass segregation for this region. 


3.2.1 Radial mass functions, Mmf 


In Fig. 5(a) the cumulative distribution of radial positions 


of the 10 most massive stars is shown by the red dashed line, 
whereas the cumulative distribution of the radial positions 
for all stars is shown by the solid black line. Due to the 
central concentration of the most massive stars, the KS test 
returns a p-value of 2 x 10““* that the two subsets share the 
same parent distribution. 

This result demonstrates that if we have confidence in 
the definition of the centre of a region, a strong mass seg¬ 
regation signature may still be seen in substructured distri¬ 
butions using the radial mass function technique. 


3.2.3 Local density ratio, Eldr 

Interestingly, the Eldr ratio does not reflect the central con¬ 
centration of the 10 most massive stars. Eldr = 0.58 due to 
the most centrally located stars being in areas of relatively 
low surface density, although a KS test returns a p-value 
of 0.26 that the massive stars have a different parent dis¬ 
tribution to the full distribution (i.e. this difference is not 
significant). That said, most people would conclude simply 
from eye that the distribution shown in Fig.|4]is mass segre¬ 
gated, even though the massive stars have low local surface 
density. In 20 realisations of this distribution, Eldr does 
not detect mass segregation in 11, and in a further 5 it finds 
inverse mass segregation. 


3 . 2.4 Group segregation ratio, Dgsr 

The overall spatial distribution has not changed between 
Figs. [ 1 ] and |4] and so the determination of dbreak and the 
subsequent identification of groups is identical to that in Sec¬ 
tion o We show the groups defined by dbreak in Fig. |6(a)[ 
noting the change of location of the 10 most massive stars 
in the full distribution (the blue open circles). The locations 
of the most massive star in each group (shown by the red 
triangles) have also changed in some cases. 

Again, we note that two of the 10 most massive stars 
from the full distribution are now no longer part of a 
group with N > 2, and instead are in a pair (located at 
{0.4, 0.1} pc). We show the mass ratio versus offset ratio for 
all groups with N > 2 stars in Fig. |6(b)| Again, the five 
groups with N ^ 10 stars are shown by the red points. 
This time only two of these five groups are mass segregated 
(Hgsr = 0.40), even though the global distribution of mas¬ 
sive stars is very mass segregated. 

In this situation - where the most massive stars in a 
region are ‘centrally’ concentrated, it is not clear what Hgsr 
is measuring. The majority of the measurements of mass 
ratio versus offset ratio do not include any of the ten most 
massive stars in the region. 


3.3 Massive stars in areas of high density 

We now change the positions of the massive stars once more, 
and swap them with stars with the highest local surface den¬ 
sities, as shown by the red points in Fig. [3 In this case we 
are shifting from a type of mass segregation related to ‘clas¬ 
sical’ mass segregation to one in which the massive stars 
are associated with the highest density regions (at least sur¬ 
face density; the volume densities of these locales would be 
unknown to the hypothetical observer). 


3.2.2 Mass segregation ratio, Amsr 


3.3.1 Radial mass functions, Mmf 


We show the Amsr ratio in Fig. 5(b) The ten most mas¬ 
sive stars have Amsr = 5.315 ®, which is significantly above 
unity. In 20 realisations, only 1 star-forming region displays 
a Amsr ratio that is not significantly above unity, with val¬ 
ues ranging from Amsr = S.OIq 5 to Amsr = 11.ll}!- 

This is completely unsurprising as Amsr is designed to 
measure exactly this type of mass segregation - the most 
massive stars being much closer to one-another than a ran¬ 
dom sample of stars would be. 


The radial mass function for the 10 most massive stars (the 
red dashed line) and the whole distribution (the black solid 
line) is shown in Fig. 8(a) As for the centrally concentrated 
cluster, the most massive stars are closer to the centre (none 
are outside of 3.3 pc) than the average stars, and a KS test 
returns a p-value of 7 x 10“® that they share the same un¬ 
derlying parent distribution. However, in 13 of our 20 reali¬ 
sations in which we draw different masses and positions for 
the stars each time, the KS test returns a p-value in excess 
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(a) A^mf (b) Amsr 


Figure 5. Three separate measures of mass segregation for the stellar distribution shown in Fig. [4| where the most massive stars are 
centrally concentrated. In panel (a) we show the cumulative distribution of the distance from the centre for the ten most massive stars 
(the red dashed line) and the cumulative distribution for all stars (the solid black line) - the mass function comparison, In panel 

(b) we show the Amsr mass segregation ratio as a function of the Nmst stars used in the subset (the lowest mass star, m^, for various 
A^mst values is shown along the top horizontal axis). Amsr = 1 no preferred spatial distribution) is shown by the solid horizontal 
red dashed line. In panel (c) we show local stellar surface density versus stellar mass (the S — m plot). The median stellar surface density 
for the ten most massive stars is shown by the righthand solid (red) horizontal line and the median surface density for all of the stars is 
shown by the blue horizontal dashed line. 



x/pc 


offset ratio 


(a) Groups identified by MST dbreak 


(b) Mass ratio versus position offset ratio for groups 


Figure 6. Mass segregation as defined by the f^csR niethod for the stellar distribution shown in Fig.[4|where the most massive stars are 
centrally concentrated. The stellar groups are identified as shown in Fig. |3(a )]and Fig. |3(by| In panel (a) we show the groups identified 
using MST dbreak- Ir groups in which there are 3 or more stars the most massive star in the group is shown by the solid red triangle. 
The positions of the ten most massive stars in the full distribution are shown by the large open blue circles (these correspond to the 
filled red circles in Fig. 0. In panel (d) we show the mass ratio of the most massive star in each group to the group median mass versus 
the ratio of the position of the most massive star to the median position of the group. Groups with ten or more stars are shown by the 
red asterisks. 


of 0.1, suggesting that the massive stars are not closer to 
the centre than the average members. This is not entirely 
surprising, as the positions in the region with the highest 
surface density may not be co-located, as is the case for one 
of the massive stars in Fig. [T] 


3.3.2 Mass segregation ratio, Amsr 

We show the measurement of Amsr for this distribution in 
Fig. 8(b) The 10 most massive stars have Amsr = 2.7lg g, 
i.e. mass segregation is present according to this measure, 
but is not as strong as in the case where we placed the most 
massive stars at the centre of the distribution. 

In many ways this is not surprising. Stars with the high¬ 
est surface density are reasonably likely to be close to one- 
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Figure 7. As Fig. ^ a fractal distribution of stellar masses ran¬ 
domly drawn from an initial mass function. However, in this case 
we have swapped the locations of the ten most massive stars 
(shown by the larger red points) with the ten stars with the high¬ 
est local stellar surface densities (as defined by Eq. 0. 



x/pc 


another (as the surface density is high), and so this artifical 
set-up will often place several massive stars close to one- 
another (here at {—0.5,—0.7} pc). This is not always the 
case, from 20 realisations of this distribution, in 11 the mea¬ 
sured Amsr was greater than unity, but less than two (and 
only marginally significant, error bars typically being around 
±0.5). 


3.3.3 Local density ratio, Eldr 


When we compare the surface density of the most mas¬ 
sive stars to the average surface density, unsurprisingly 
the most massive stars have much higher median values. 


as shown in Fig. 8(c) Here, the most massive stars have 
E = 64.1 starspc“'^, compared to E = 13.1 starspc“^ for 
the region average. Eldr = 4.9, and a KS test returns a p- 
value of 8 X 10“^ that they share the same underlying parent 
distribution. 

This is exactly as expected as the set-up is such that 
Eldr should find mass segregation by its definition of it. 


3 . 3.4 Group segregation ratio, Hgsr 

As in Section S/Ud shOAVn in Fig. C^break 1® th© 

same as that calculated in Section ED because the spatial 
distribution has not changed. The groups identified by MST 
dbreak are shown in Fig. 9(a)[ and the most massive star in 
each group is shown by the red triangle. The 10 most mas¬ 
sive stars in the distribution are shown by the blue circles. 
This time, none of these 10 massive stars are not in groups, 
but we have a significant problem that one group contains 9 
of them (and so 8 will be discarded from the analysis). When 
we determ ine whether that grou p is mass segregated accord¬ 
ing to the lKirk fc Mvei^ (l201lll method, we are effectively 
ignoring the positions of 8 of these stars. This time, three 
of the five groups with V ^ 10 are mass segregated, and 


the largest group (containing 9 of our most massive stars 
in the full distribution) is also mass segregated according to 
Amsr and Eldr- However, the other large group (containing 
the single massive star) i s not mass segregated according to 
the lKirk fc Mve^ 1I2OIII) method, but it is with Amsr and 
Eldr. For the groups with N > 2, Hqsr = 0.59, and for 
N ^ 10 Hgsr = 0.60, i.e. more groups than not are mass 
segregated according to this method. 


4 DISCUSSION 


When attempting to find ‘mass segregation’ in a region it is 
absolutely critical to clearly define what is meant by ‘mass 
segregation’. Confusion between apparently contradictory 
results for ‘mass segregation’ between different methods oc¬ 
curs because the differe nt methods are searching for d iffer- 
ent things. For example, iMaschberger fc Clarkel (I 201 J) find 
mass segregation according to El dr in the hydrodyna mical 
simulations of star formation from iBonnell et ;rD ll2008h . but 
do not find mass segregation with Amsr. We contend that 
this is not condradictor y, rather just differ ent definitions of 
‘mass segregation’ (e.g. Parker et al.ll2014h . 

A definition of mass segregation based on relaxation and 
equipartition in a dynamically old system is one in which 
the most massive stars are closer together than expected by 
random chance. It is this definition that is proped by radial 
mass function methods (A4mf) and Amsr- In searching for 
this type of mass segregation Amsr is more useful as it does 
not require a centre to be defined and can deal with complex 
(substructured) distributions. 

The Eldr method defines mass segregation differently 
- in this case mass segregation is that the most massive 
stars are preferentially in regions of higher surface density 
than random. Whilst this method does not measure ‘mass 
segregation’ in the classical sense, it is extremely useful for 
probing the past dynamical history of a star-forming re¬ 
gion, as the most massive stars sweep up retinues of low- 
mass stars during the two -body relaxat i on of i nitially dense 
(>100M(7) pc~^) r egions (IParker et al.l I20l4 IParked I20l4 
I Wright et al.ll20lj ). 

One way of avoiding confusion between Amsr and Eldr 
is to make the definition of mass segregation more stringent, 
for example that the most massive stars should be globally 
more concentrated and be at the centre of individual groups. 
However, this requires the somewhat arbitrary definition of 
groups, which is arguably impossible if all the stars formed 
in the same star formation ‘event’ in the same molecular 
cloud, and so any boundary between groups is necessarily 
artificial. Furthermore, as we have seen in Section [3.21 a spa¬ 
tial distribution that few would argue is not mass segregated 
would fail this definition. 

In this context it is unclear to the authors what ex¬ 
actly Hgsr is searching for, or what definition of ‘mass seg¬ 
regation’ it involves. We have also identified a number of 
problems with the Hgsr method which we feel makes it un¬ 
suitable for finding ‘mass segregation’. 

Firstly, it is unclear if the group identification is in any 
way finding ‘real’ groups. 

Secondly, the identification of groups discards any stars 
that are not in an N > 2 group (or higher N, depending 
on the number of stars per group as defined in a particular 
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(a) A^mf 


(b) Amsr 



Figure 8. Three separate measures of mass segregation for the stellar distribution shown in Fig. [7| where the most massive stars are 
in the areas of highest stellar surface density. In panel (a) we show the cumulative distribution of the distance from the centre for the 
ten most massive stars (the red dashed line) and the cumulative distribution for all stars (the solid black line) - the mass function 
comparison, AImf- In panel (b) we show the Amsr mass segregation ratio as a function of the A’mst stars used in the subset (the 
lowest mass star, mi,, for various A’mst values is shown along the top horizontal axis). Amsr = 1 (i-®. no preferred spatial distribution) 
is shown by the solid horizontal red dashed line. In panel (c) we show local stellar surface density versus stellar mass (the S — m plot). 
The median stellar surface density for the ten most massive stars is shown by the righthand solid (red) horizontal line and the median 
surface density for all of the stars is shown by the blue horizontal dashed line. 
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(a) Groups identified by MST dbreak 


(b) Mass ratio versus position offset ratio for groups 



Figure 9. Mass segregation as defined by the Hgsr method for the stellar distribution shown in Fig. [T] where the most massive stars 
are in areas of highest stellar density. The stellar groups are identified as shown in Fig. |3(a)| and Fig. |3(b)| In panel (a) we show the 
groups identified using MST dbreak- Ir tbe groups in which there are 3 or more stars the most massive star in the group is shown by 
the solid red triangle. The positions of the ten most massive stars in the full distribution are shown by the large open blue circles (these 
correspond to the filled red circles in Fig. 0. In panel (d) we show the mass ratio of the most massive star in each group to the group 
median mass versus the ratio of the position of the most massive star to the median position of the group. Groups with ten or more stars 
are shown by the red asterisks. 


analysis. iKirk fc Mver j[201ll : lKirk et al.ll2014l l. This ignores 
many stars, removing them from further analysis, even if 
they are amoung the most massive stars in the region. 

Thirdly, once groups have been identified the method 
only considers the most massive star in that group, discard¬ 
ing information on the masses of any other stars. 

Forthly, groups are assumed to have a ‘centre’ from 
which distances can be measured, essentially performing a 


‘radial mass function’ approach based on a single massive 
star in a small-V subset of the total population. 

Finally, it is worth noting that small-V statistics can 
play an important role in obtaining any information from 
the flcsR method. In the examples we showed above there 
are only five groups with N > 10. For no ‘mass segregation’ 
we would expect 2 or 3 to show no signal, however it would 


© 2015 RAS, MNRAS OOO.fTlfTSl 
















12 R. J. Parker & S. P. Goodwin 


not be unusual for only 0 or 1 to show no signal, or 4 or 5 
to. 

Variations on this final point is important for all meth¬ 
ods. A positive signal for ‘mass segregation’ in-and-of-itself 
may not tell us much. As we saw in the example random dis¬ 
tribution we used above (Section 13. 1.211 . Amsr found mass 
segregation at 2-cr significance. This is a result we would 
expect l-in-20 times, and examining our ensemble of simu¬ 
lations we Hnd this is indeed the case (and some show ‘in¬ 
verse mass segregation’ in which Amsr < !)• This ‘random 
noise’ effect has b een seen in ensembles of simulations (see 
iParker et al.ll2015l l. 

Based on this, it is quite possible that small signatures 
of ‘mass segregation’ such as the appare ntly inverse mass 
segregation found by IParker et al.l (120111 1 in Taurus might 
well have been over-interpreted and are quite possibly con¬ 
sistent with a random distribution of the most massive stars. 
This highlights the requirement for more than one technique 
to be applied to any search for mass segregation in an ob¬ 
served region. 


5 CONCLUSIONS 


We have experimented with four methods used to find ‘mass 
segregation’ : the r adial mass functi on metho d Mmf (e.g. 


Saear et al.ll 19881: Sabbi et al.l 2008h. Amsu (I Allison et al.l 

2009), Eldr 

Maschbereer & Clarkel boilll. and Hnsu 

(Kirk & MversI 

201lf). Our results can be summarised as 


follows. 


(i) Only in smooth, spherical, centrally concentrated 
distributions (e.g. Plummer spheres) do all methods find 
‘mass segregation’. In more complex, substructured distri¬ 
butions different methods can find different things because 
they define ‘mass segregation’ differently. 

(ii) Only Amsr measures ‘classical’ mass segregation 
where the massive stars are concentrated in particular re¬ 
gions without having to define a cluster centre. 

(hi) The radial mass function method A4mf searches 
for ‘classical’ mass segregation, but requires a centre to be 
defined and then assumes spherical symmetry. 

(iv) Eldr measures a different ‘mass segregation’ where 
the massive stars are in regions of higher than average sur¬ 
face density without having to define a cluster centre. The 
massive stars may, or may not, also be concentrated to¬ 
gether. 

(v) flcsR finds groups that may, or may not, be 
physically important and then defines a ‘centre’. In doing 
so it can exclude very significant information on some of 
the most massive stars in a region (sometimes excluding 
them from the analysis entirely). 

We conclude that of the methods currently in use, by 
far the most useful are Amsr and Eldr. They use all of the 
information on all of the stars in a region without assum¬ 
ing anything about the spatial distributions. We reiterate, 
however, that they measure different definitions of ‘mass seg¬ 
regation’ and so should be used in tandem. 

Finally, we note that marginal signals of mass segrega¬ 
tion (as found by any method) in observed star-forming re¬ 
gions may not have anything to do with the physics of star 


formation, and any analysis should be accompanied by a 
suite of simulations of synthetic regions like those presented 
here. In a future paper, we will also examine the poten¬ 
tially significant and serious problems of analysing projected 
distributions and attempting to extract information on the 
three dimensional properties. 
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